Systematic tight-binding analysis of ARPES spectra of transition-metal oxides 
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We have performed systematic tight-binding (TB) analyses of the angle-resolved photoemission 
spectroscopy (ARPES) spectra of transition-metal (TM) oxides AMO3 (M = Ti, V, Mn, and Fe) 
with the perovskite-type structure and compared the obtained parameters with those obtained from 
configuration-interaction (CI) cluster-model analyses of photoemission spectra. The values of e^ — e p 
from ARPES are found to be similar to the charge-transfer energy A from O 2p orbitals to empty TM 
3d orbitals and much larger than A — U/2 (U: on-site Coulomb energy) expected for Mott-Hubbard- 
type compounds including SrVOs. e<i — e p values from ab initio band-structure calculations show 
similar behaviors to those from ARPES. The values of the p — d transfer integrals to describe the 
global electronic structure are found to be similar in all the estimates, whereas additional narrowing 
beyond the TB description occurs in the ARPES spectra of the d band. 

PACS numbers: 79.60.-i, 71.15.Mb, 71.28.+d, 71.30.+h 



3d transition-metal (TM) oxides have attracted a lot of 
interest in these decades because of their intriguing phys- 
ical properties such as metal-insulator transition, colos- 
sal magnetoresistance, and the ordering of spin, charge, 
and orbitals [l[ . Photoemission spectroscopy has greatly 
contributed to the understanding of the electronic struc- 
tures of these materials. The local electronic structure of 
strongly correlated 3d states of the TM atoms hybridizing 
with ligand p orbitals has been understood in the frame- 
work of the Zaanen-Sawatzky-Allen (ZSA) diagram [2( 
based on the configuration-interaction (CI) cluster-model 
analyses of core-level and valence-band angle-integrated 
photoemission spectra 0, 0, @ • On the other hand, how 
the energy band structure of the periodic lattice emerges 
from the correlated local electronic structure has been an 
extremely difficult problem as has been extensively dis- 
cussed in the case of /-electron systems [1, 0], and no 
clear scenario has been established so far. On the exper- 
imental side, rather surprisingly Sarma et at Q have 
found that the angle-integrated valence-band photoemis- 
sion spectra of various perovskite-type 3d TM oxides can 
be well reproduced by band-structure calculation based 
on the local spin-density approximation in spite of the 
strong electron correlation. 

In order to address the issue of a proper description 
of the periodic systems with strong correlation, the ex- 
perimental determination of band dispersions by angle- 
resolved photoemission spectroscopy (ARPES) has been 
highly required, but such data have been largely lim- 
ited to low dimensional compounds 0, [l(J E3 of limited 
chemical variety. For three-dimensional perovskite ox- 



ides, which have a wide variety of chemical compositions, 
even when bulk single crystals are available, they do not 
have cleavage planes required for ARPES measurements. 
Recently, technological development has made it possible 
to grow high-quality single-crystal thin films using the 
pulsed laser deposition method, and a setup has been 
developed for their film growth followed by in- situ pho- 
toemission measurements [Til. K$L Hi . 15 1. Since ARPES 
data have accumulated for various 3d TM oxides of per- 
ovskite structures [SrTi0 3 (STO) J3, SrV0 3 (SVO) 
0, Lai-^Sr.MnOa (LSMO) 0, Qjf PiwCa^MnOa 
(PCMO) [H and Lai_ x Sr x Fe0 3 (LSFO) 0] and have 
been analyzed using a parameterized nearest-neighbor 
tight-binding (TB) model [H [H, 0, HH, we are in a 
position to look at the systematic variation of the elec- 
tronic structure of the periodic lattice of 3d TM atoms 
and derive empirical rules, as has been done for the lo- 
cal electronic structure by ZSA In particular, since 
there is no established theoretical description for corre- 
lated periodic systems, it is highly desirable to obtain 
relationship between parameters in CI theory and those 
describing ARPES band structures. In this paper, we 
examine the TB parameters that describe the ARPES 
band dispersion in order to deduce systematic variations 
of the energy level (e^) of TM 3d orbitals relative to 
that (e p ) of O 2p orbitals and the p — d transfer inte- 
grals [Slater-Koster parameter (pda)] and compare them 
with those derived from the CI cluster-model analyses of 
angle-integrated photoemission spectra, namely, {pda), 
the charge-transfer (CT) energy from the O 2p orbitals 
to the empty TM 3d orbitals denoted by A, and the aver- 



2 



age 3d — 3d on-site Coulomb interaction energy denoted 
by U . For comparison, we have also fitted the ab initio 
band structure to the same TB model, as has been done 
in Ref. (26[. The values of e d — e p from ARPES were 
found to be similar to A rather than A — U/2 expected 
for Mott-Hubbard (MH) type compounds. The values 
of |(p<i<7)| (determined from the O 2p band dispersions) 
were found to be similar in these estimates, but addi- 
tional narrowing in the d band (by a factor of 0.3 — 0.5) 
was seen in ARPES beyond the TB description. 

All the ARPES data analyzed in this paper were mea- 
sured using a photoemission spectroscopy system com- 
bined with a laser molecular beam epitaxy chamber at 
beamline BL-1C of the Photon Factory, KEK [T3|. De- 
tails of the experimental conditions for the respective 

materials were described in Refs. [H, E3, 0, H, EE EH . 

The ab initio band-structure calculations were performed 
using the linearized augmented plane wave method im- 
plemented in the WIEN2K package [13] • The exchange 
and correlation effects are treated within the generalized 
gradient approximation [28| • We assumed a cubic param- 
agnetic (PM) state, and the lattice constants were taken 
as 3.91A for STO, 3.84A for SVO, 4.04A for LaMn0 3 
(LMO), and 4.02A for LaFe0 3 (LFO). We fitted the 
ARPES and ab initio results to parameterized nearest- 
neighbor TB band structures [3 HE HI- We as- 
sumed the PM state for the ARPES results of STO and 
SVO, the ferromagnetic (FM) state for LSMO, the C- 
type antiferromagnetic (AF) state for PCMO, and the 
G-type AF state for LSFO. The size of the matrix for 
the calculations were 14 x 14 (TM 3d orbitals: 5, oxygen 
2p orbitals: 3 x 3 = 9) in the cases of PM and FM states, 
and 28 x 28 in the cases of AF states. The effects of fer- 
romagnetism and antiferromagnetism were treated phe- 
nomenologically by assuming an energy difference AE 
between the spin-up and spin-down TM sites. Parame- 
ters to be adjusted are AE, the energy level of the TM 3d 
orbitals relative to that of the O 2p orbitals, ed — e p , and 
Slater-Koster parameters (pda), (pdir), (ppcr), and (ppn) 
[iij ]. In the cases of the FM and AF states, ed is the av- 
erage of the up-spin d level ed — AE/2 and the down-spin 
d level ed + AE/2. The details of the matrix elements for 
the calculations were described in Refs. [23j,|25j . 

In the CI cluster-model calculations [J |j7 !, H3] , a 
TM06 octahedral cluster is considered. In this model, 
the ground state is described as 

* 9 = ai \d n ) + a 2 \d n+1 L) + a 3 \d n+2 L 2 ) + ■■■. 

The final state for the emission of a TM 2p core electron 
is given by 

*/ = Pi\cd n } + I32\cd n+1 L) + (3 3 \cd n+2 L 2 ) + • • • , 

where c denotes a TM 2p core hole. The final state for 
the emission of a TM 3c? electron is given by 

*/ = 7iK _1 ) + l2\d n L) + -f 3 \d n+1 L 2 ) + ■■■. 
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FIG. 1: (Color online) Tight-binding fit to the ARPES spec- 
tra [TtJ (a) and ab initio band-structure calculation (b) for 
SrV03. In panel (a), bright parts correspond to energy bands. 



In this work, we have taken the parameter values of STO 
and SVO from Ref. '30], and those for LSMO and LSFO 
from Refs. 
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32j, where these values have been ob- 
tained by fitting the TM 2p core-level and valence-band 
photoemission spectra. 

In Fig. [TJ we show an example of the TB fit to the 
ARPES spectra and ab initio band-structure calculation. 
Figure [TJ (a) shows a fit to the ARPES spectra of SVO 
[171 ]. The bands below ~ —3 eV are mainly from O 2p 
states, and those between ~ —2 eV and Ep are mainly 
from V 3g! states. The dispersions of the O 2p bands were 
well reproduced by the calculation, but the widths of the 
V 3d bands were about 50% of the calculated ones. Such 
a mass renormalization is also observed in the ARPES 
spectra of LSMO 18 1 and high-Tc cuprates Figure 
UJ (b) shows a TB fit to the ab initio band-structure cal- 
culation of SVO. The parameters ed — e p — 3.6 eV and 
(pda) = —2.2 eV give a fairly good fit to the ab initio 
calculation. 

Figure [2] shows comparison of the variation of | (pda) 
[panel (a)] and — e p [panel (b)] for ARPES and ab 
initio compared with parameter values deduced from CI 
cluster-model analysis. As seen from Fig. [2] (a), the val- 
ues of (pda) are similar in all the three estimates. This 
means that, as far as the strength of p — d hybridization 
is concerned, the values from the local CI cluster-model 
analysis can be used to describe the band structures of 
periodic systems. The behavior of — e p from ARPES 
follows A from CI analysis, as shown in Fig. [2] (b). The 
£<i — £p values from ab initio calculations are similar to or 
a little (by ~ 1 eV) smaller than those from ARPES. This 
reflects the well-known tendency of ab initio calculations 
to underestimate the band gaps and the binding energies 
of the O 2p bands. One typical example is the underes- 
timation of the values of the bad gaps in STO [33j. The 
indirect (R — > T) and direct (r — > T) gaps are calculated 
to be 1.9 eV and 2.2 eV, while the experimental values 
are 3.25 eV and 3.75 eV. However, as this difference is 
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FIG. 2: (Color online) Material dependences of the p—d trans- 
fer integral (\(pda)\) (a) and the energy level of TM d relative 
to O 2p (ed — £p) (b) deduced from ARPES and ab initio calcu- 
lations compared with those from CI calculation. Values from 
ab initio calculations in Ref. |26| has been also plotted. The 
CI parameter values of STO and SVO are taken from Ref. [30| . 
and those for LSMO and LSFO are from Refs. [USl. 
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rather small, one may conclude that ab initio calculations 
can describe the angle-integrated photoemission spectra 
of these materials fairly well [8|. It is remarkable to see 
that, except for the common shift of the O 2p band, the 
experimental ed — e p of complicated origin is well repro- 
duced by the simple ab initio calculation. 

For the d° compound (STO), the value of ed — e p is 
expected to be given by A. In the case of MH-type com- 
pounds such as SVO, ed — e p is expected to be given by 
A — U/2 [201 . However, Fig. [5] shows that the ed — e p 
values from ARPES are closer to A than A — U/2 for all 
compounds including SVO. One should note here that 
the values of A shows strong dependence on TM atoms 
and changes from ~ to ~ 6 eV, whereas those of U do 
not strongly depend on TM atoms and are as large as 
5-8 eV liEi 

To understand the difference in ed — e p between ARPES 
and the simple MH picture even for SVO, we consider the 
single-particle energy diagrams for MH type compounds, 
where A > U, [panel (a)] and CT type compounds, where 
A < U, [panel (b)] and the assignments of band charac- 
ters as shown in Fig. [3] In a MH-type compound, the 
TM 3d states are located above the O 2p band, and the 
p—d hybridization simply causes upward and downward 
shifts of the p and d levels, respectively, leading to little 
difference between ARPES and CI. Therefore, in pure 
MH-type compounds, the relationship ed — e p = A — U/2 
holds, which is inconsistent with our results. In a CT- 
type compound, the effect of p — d hybridization is com- 
plicated, creating new states in the occupied part, i.e., in 
the electron removal spectrum. That is, a charge-transfer 
satellite of eP -1 final-state character appears below the 
O 2p band and the split-off state of the Zhang-Rice (ZR) 
states (d n L final-state character) are formed above the 
O 2p band as a result of p — d hybridization, as shown 
in Fig. [3] (b). In the case of a d 9 system (Cu 2+ ), the ZR 



FIG. 3: (Color online) Energy diagrams for Mott-Hubbard- 
type (a) and charge-transfer- type (b) compounds. Sat: 
charge-transfer satellite (d n ); ZR: Zhang- Rice singlet-like 
states (d n L). Here, ed — e p is the 3d level relative to O 2p, A 
is the charge-transfer energy from O 2p to the empty TM 3d, 
and U is the 3d — 3d on-site Coulomb interaction energy. 



state becomes a spin singlet, well known as the ZR singlet 
[3~ij ] . The emergence of the charge-transfer satellite and 
the ZR state can be understood from the Anderson im- 
purity model, where the discrete TM 3d impurity states 
interact with the O 2p continuum states through strong 
p — d hybridization and two new split-off states are cre- 
ated. When we fit the ARPES results, we consider the 
ZR state as "effective" TM 3d bands. This explains why 
the ed — e p values in ARPES becomes considerably larger 
than A — U/2 in the CT-type compounds as shown in 
Fig. OJb). The difference of these two is large also in 
SVO, which has long been considered as a typical MH 
compound. This result shows that SVO is not a pure 
MH compound and hybridization between O 2p and V 
3ci states is comparable to CT-ty pe compounds, as al- 
ready suggested in Fig. 15 in Ref . 13011 and concluded by 
the recent CI calculation in Ref. [35j. 

The similarity of the ed — e p values in ab initio and 
ARPES shows that in ab initio ed — e p is similar to A not 
to A — U/2. This result means that the d bands in ab 
initio calculations are not the pure d bands before p — d 
hybridization but the ZR state after being considerably 
hybridized with O 2p states. 

Finally, we would like to comment the limitation of our 
analyses of the ARPES spectra. The values of |(pder)|, 
which describes the global electronic structure, were sim- 
ilar in the cases of ARPES, ab initio, and CI, but the 
ri-band narrowing observed in ARPES were about 1 /2 or 
1/3, which cannot be described in the TB band-structure 
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calculation. This is the effect of electron correlation at 
the TM sites. 

In conclusion, we have compared the parameters ob- 
tained from TB fit to the ARPES spectra of TM oxides 
AM0 3 (M = Ti, V, Mn, and Fe) and those obtained by 
TB fit to the ab initio band-structure calculation with 
those from CI cluster-model calculations. The — e p 
values obtained by TB fit to ab initio band-structure cal- 
culations were a little smaller than those in ARPES, re- 
flecting the general tendency of ab initio calculations to 
underestimate the binding energies of O 2p bands. The 
£d — £p values from ARPES were significantly larger than 
A — U/2 for all the compunds including SVO, which has 
been considered as a typical MH system. We attribute 
this trend to the ZR-band origin of the effective d band. 
The values of |(pd<r)| were found to be similar in all the 
three estimates, showing the validity of using parame- 
ters obtained from the local cluster or impurity Ander- 
son models to describe the global electronic band struc- 
ture of the TM oxides. However, the narrowing of the 
("effective") d band seen in ARPES was beyond the TB 
description. 
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